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Abstract 



In this thesis, we perform the lattice QCD analysis via the energy-momentum 
component of gluons. By introducing the momentum cutoff to the link variable, we 
investigate which energy-momentum components of gluons induce each QCD phe- 
nomenon. We use the Landau gauge for the most part of the lattice QCD analysis. 
In lattice QCD, we analyze color confinement, spontaneous chiral symmetry break- 
ing, topological charge, and the related topics. We also discuss several comparisons 
with effective theories. As for color confinement, wc calculate the quark- ant iquark 
potential, the color flux tube, and meson masses. From quantitative analysis, we 
flnd that color confinement is induced by the low-momentum component below 1.5 
GeV. As for spontaneous chiral symmetry breaking, we calculate the chiral conden- 
sate and the Dirac spectrum. Spontaneous chiral symmetry breaking is induced by 
the broad low-momentum component which ranges even above 1.5 GeV. The present 
result suggests that color confinement and spontaneous chiral symmetry breaking 
are induced by somehow different energy-momentum components of gluons. As for 
topological charge, we calculate the topological charge density and the Dirac zero 
mode. Topological charge is induced by the broad energy-momentum component, 
which is similar to spontaneous chiral symmetry breaking. 
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Chapter 1 
Introduction 



Physics and scale 

In physics, the scale is ranged from the smallest size (the Planck scale) to the 
largest size (the universe) . There exist a huge number of physical phenomena. These 
phenomena are described in various ways depending on their scales. The small-scale 
phenomena are described by elementary particle physics, nuclear physics, and so on, 
and the large-scale phenomena are described by cosmology, astrophysics, and so on. 
The scale is one of the most fundamental and important concepts in physics. 

In quantum chromodynamics (QCD), there is a characteristic property describing 
the nontrivial appearance of a scale. The QCD Lagrangian is composed of gluons 
and quarks as 

£qcd = -\i:t[F^''F,,] + g[z7M^^ - m]q . (1.1) 

This Lagrangian has no dimensional parameter except for quark masses. Thus, 
classical QCD is scale invariant in the chiral limit m — )■ 0. After the quantization, 
however, scale invariance is violated by the trace anomaly, and a scale appears in 
a nontrivial manner. This is called as dimensional transmutation. In this way, 
many dimensional quantities are created in real QCD, such as "mass gap" , vacuum 
condensates, and masses of hadrons. 

As a result of the dimensional transmutation, the QCD running coupling constant 
is a function of the energy scale (Fig. II. ip . At the one- loop perturbation theory, the 
QCD running coupling constant is 

4vr -4vr/5oln(QVA^cD)' 

where /3o = (HiVc - 2Nf)/A8TT^ [HE]. For this reason, the behaviors of QCD phe- 
nomena are completely different among different energy scales. At high energy or 
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Figure 1.1: QCD running coupling constant as(Q) [3j. 



short distance, perturbative QCD is valid due to the asymptotic freedom. At low en- 
ergy or long distance, the interaction is strong coupling and nonperturbative effects 
are important. The appearance of the energy scale enriches QCD phenomenology. 

In this study, our question is "which energy scale is relevant for each QCD phe- 
nomenon?" It is often said that Aqcd is the typical scale of QCD. However, our aim 
is not to determine Aqcd- Aqcd is not a quantitative scale of each phenomenon, 
but a typical scale of the whole theory. Actually, even above Aqcd, the running 
coupling constant is still large and the nonperturbative effect is important. Since 
the energy scale is abstract concept, we have to define the energy scale in a more 
concrete way. We define the relevant energy scale for a QCD phenomenon as the 
energy-momentum component of gluons inducing the phenomenon. Thus, our ques- 
tion is rephrased as "which energy-momentum components of the gluon field induce 
each QCD phenomenon?" By definition, such a relevant energy scale is different 
from Aqcd- It is nontrivial whether the relevant energy scales are the same or dif- 
ferent between different phenomena, e.g., color confinement and spontaneous chiral 
symmetry breaking. For a scale-dependent theory like QCD, the knowledge of the 
relevant energy scale would be useful for understanding physical phenomena. 

To answer this question, we employ the numerical simulation of lattice QCD. 
Lattice QCD is one of the most promising calculations for nonperturbative QCD. 
In this thesis, we introduce a lattice QCD framework to determine which energy- 
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momentum components of the gluon field induce QCD phenomena. Using the in- 
troduced framework, we discuss the relevant energy scales, systematically, quanti- 
tatively, and nonperturbatively. 

Color confinement and spontaneous chiral symmetry breaking 

We focus on two of the most significant phenomena in QCD; color confinement 
and spontaneous chiral symmetry breaking. 

In hadron phase, a color source, such as a quark or an antiquark, is confined 
in hadrons. Color confinement symbolically reflects the strong interacting nature 
of QCD. Phenomenologically, it has long been known that color confinement is 
well described by a "string" or "flux tube" |4j. The analytical derivation of color 
confinement has not yet been achieved because of the nonperturbative and non- 
Abelian nature of QCD. It is not only an important theme in physics but also 
an unsolved problem in mathematics. Lattice QCD studies have shown that color 
confinement is indeed realized, i.e., the Wilson loop obeys the area law [5l[6]. Color 
confinement is one of the main subjects in lattice QCD [3 |8]. 

Although the QCD Lagrangian possesses chiral symmetry except for the quark 
mass term, this symmetry is spontaneously broken into its subgroup as 

S\J{Nf)L X S\J{Nf)R S\J{Nf)v. (1.3) 

Spontaneous chiral symmetry breaking is one dominant origin of mass in our world, 
and it is remarkably important in various aspects of hadron physics [9l [TOl [HI [T2] . 
Chiral symmetry itself is symmetry of quarks, not of gluons. However, spontaneous 
breaking is dynamically induced by the nonperturbative interaction of gluons. The 
gluon dynamics is essential for spontaneous chiral symmetry breaking. The relation 
between the eigenmode of quarks and chiral symmetry breaking is known as the 
Banks-Casher relation [13]. On the other hand, the relation between the energy- 
momentum component of gluons and chiral symmetry breaking is nontrivial. 

The connection between color confinement and spontaneous chiral symmetry 
breaking is often discussed in QCD [HI [TSl HSl [IZ]. At finite temperature, de- 
confinement and chiral symmetry restoring phase transitions occur at almost the 
same temperature [T8| \T9\ 120] . We explore this connection in the viewpoint of the 
energy-momentum components of the gluon field at zero temperature. 

Outline of the thesis 

The outline of this thesis is as follows. In Chapter 2, we briefly review basics 
of lattice QCD, and introduce our framework to analyze the relation between QCD 
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phenomena and the energy-momentum component of the gluon field. Using this 
framework, we determine the relevant energy scale in SU(3)c lattice QCD. To inves- 
tigate color confinement, we analyze the static quark-antiquark potential in Chapter 
3, the color fiux tube in Chapter 4, and the meson masses in Chapter 5. To inves- 
tigate spontaneous chiral symmetry breaking, we analyze the chiral condensate in 
Chapter 6. and the Dirac spectrum in Chapter 7. Each simulation setup is shown 
in each chapter. Finally, Chapter 8 is devoted to summary. 
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Chapter 2 



Lattice QCD Formalism 



2.1 Basics of lattice QCD 



Partition function 

Lattice QCD is one of the most powerful tools to analyze nonperturbative QCD 
phenomena |2T1I221E3]- In lattice QCD, the space-time is discretized as a hypercubic 
lattice with the lattice spacing a. In the following, we consider the four dimensional 
space-time in the Euclidean metric. 

The QCD partition function is given by 



The lattice action is constructed so as to reproduce the continuum action in con- 
tinuum limit a — )■ 0. The expectation value of a physical quantity is computed by 
Monte Carlo simulation as 



The ensemble average is taken over the gauge configurations, which are generated 
with the weight function ~ exp(— S'gauge — •S'quark)- The fermion action in the weight 
function is sometimes dropped for simplicity, which is called as the quenched ap- 
proximation. The quenched approximation corresponds to neglecting the quark loop 
effects. 




(2.1) 




(2.2) 



(2.3) 
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Gauge field 

In lattice QCD, the SU(A^c) gauge field is represented as the link variable 

U^{x) = e"'3^'^^^\ (2.4) 
The minimum gauge invariant combination of the link variable, 

U,,{x) = U,{x)U,{x + fi)Ul{x + 0)Ul{x), (2.5) 
is called as a plaquette. The simplest lattice gauge action is written by this plaquette 



as 



5'gauge = ^s{x) (2.6) 

X 

six) = /3 5^(l-^ReTrf/^,(a:)) (2.7) 

^ - 2^. (2.8) 

This lattice gauge action corresponds to the continuum gauge action in continuum 
limit a — )■ 0. 

The expectation value of a gauge variant operator is automatically zero without 
gauge fixing because of Elitzur's theorem [21]. To calculate a gauge variant operator, 
gauge fixing is necessary. Instead of the Faddeev-Popov method, gauge fixing is 
realized by imposing the gauge condition numerically. For example, the Landau 
gauge is implemented by globally maximizing the quantity 

Fi[f/] = ^^ReTr[/^(x), (2.9) 

X fl 

by the gauge transformation. In continuum limit, this condition is equivalent to 
minimizing 

I rf^xTr{A^(x)^}, (2.10) 
and it is a sufficient condition for the local condition df^Af^{x) = 0. 

Fermion field 

The fermion action is the bilinear form of the quark field, 

Vark = 5Zg(x)/^g(z/). (2.11) 

x,y 



2.2. MOMENTUM CUTOFF 
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If one naively discretizes the continuum Dirac operator D = "y^d^ — ig'-f^A^, the 
lattice Dirac operator suffers from the doubling problem [25]. Several ways have 
been proposed to avoid the doubling problem. 

One familiar lattice fermion is the Wilson fermion [26]. The Wilson Dirac oper- 
ator is 

D = J]{(l-7^)t/^(a;)5,+^,, + (l + 7^)?7^(a;-/i)5,_^,,}, (2.12) 

1 

K = . 2.13 

2ma + 8 ^ ^ 

The Wilson Dirac operator includes the O(a^) term which breaks chiral symmetry 
explicitly. The degree of freedom of flavor can be introduced without restriction. 

Another familiar lattice fermion is the staggered fermion f27\ [25] . The staggered 
Dirac operator in the spinorless basis is 

D = m6^,y + ^'^r]^{x)[Uf,{x)6^+f,^y-Ul{x- fi)6^-f,J , (2.14) 

T^^(x) = (_i)^i+-+^M-i . (2.15) 

The staggered Dirac operator describes the degenerate Nf = 4 Dirac fermion in 
continuum limit. At finite lattice spacing, however, it includes an O(a^) term which 
violates the flavor symmetry, and the symmetry is reduced toU(l)yxU(l)^. Com- 
pared to the Wilson fermion, the staggered fermion has the 11(1)^ subgroup of the 
chiral symmetry which is exact on the lattice, but is restricted in the number of 
flavor. 

The above two fermions break the full chiral symmetry. As a result, the anti- 
commutation relation {D,'-/^} = 0, which holds in continuum QCD, is not satisfied. 
There are other lattice fermions with better chiral property, such as the overlap 
fermion [29], the domain- wall fermion [30], and so on. 



2.2 Momentum cutoff 

To investigate the relevant energy scale, we introduce a lattice framework to 
remove some region of momentum space by a momentum cutoff. By observing 
how a physical quantity is affected by the momentum cutoff, we can determine the 
role of the removed energy- momentum component. This concept has been applied, 
for example, in the Swinger-Dyson approach [21]. We introduce this concept to 
lattice QCD, and analyze the energy-momentum component of the gauge field. The 
framework is applicable to both quenched and full QCD in the same way. 

The framework is formulated as the following five steps. 
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Step 1 

The SU(3)c link variable U^{x) is generated by Monte Carlo simulation. Bound- 
ary conditions for the link variable are taken to be periodic. As explained below, 
the link variable must be fixed with a certain gauge. In the following, we mainly 
use the Landau gauge for the numerical calculation, unless otherwise stated. 



Step 2 

The momentum-space link variable U^lp) is obtained by the Fourier transforma- 
tion, as 

Utiip) ^ yYl ^1^^^^ ®^P(^ ^Pv^i^), (2.16) 

X V 

where V is the lattice volume. The momentum-space lattice has the same structure 
as the coordinate-space lattice. Its lattice volume is V , and the boundary conditions 
are periodic. The momentum-space lattice spacing is given by 

= (2.17) 

where L is the number of lattice sites in each direction. The momentum-space lattice 
spacing corresponds to the minimum unit of momentum, and it has mass dimension. 
The first Brillouin zone of the momentum space is (— vr/a, vr/a]. 



Step 3 

Some components of U^ip) are removed by introducing a momentum cutoff. In 
the cut region, the momentum-space link variable is replaced by the free-field link 
variable, 

I (other region). 



where 

IX vx-/ y 

The concrete form of the momentum cutoff can be taken arbitrarily. For example. 



^'"'b) = ^ E 1 exp(^E.P.^.) = V- (2-19) 



X 



one natural choice is the cutoff by the four- momentum length yp^ = y ^^P/^P/x, 
which corresponds to a simple momentum cutoff in continuum theory. The ultravi- 
olet (UV) cutoff Auv is introduced as 

U,{p) (v^<Auv) 



' (VP > Auv), 



2.2. MOMENTUM CUTOFF 
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Figure 2.1: The schematic figures of momentum space. The shaded regions are the cut 
regions by the UV cutoff Auv (left) and the IR cutoff Air (right). The momentum-space 
lattice spacing is ap = In /La. 



and the infrared (IR) cutoff Air is introduced as 



< Air) 
^ > Air) 



(2.21) 



The schematic figures are shown in Fig. 12. 1[ 



Step 4 

The coordinate-space hnk variable with the momentum cutoff is obtained by the 
inverse Fourier transformation as 

K(^) = $^t/^(p)exp(-z5^p,x,). (2.22) 

p V 

Since U'Jyx) is not an SU(3)c matrix in general, U'^{x) must be projected onto an 
SU(3)c element U^{x). The projection is realized by maximizing the quantity 

ReTr[{f/^^(x)}tf/;(a:)]. (2.23) 

By this projection, we obtain the coordinate-space link variable f/^(a;) with the 
momentum cutoff, which is an SU(3) matrix and has the maximal overlap to U'^{x). 
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Step 5 

The expectation value of the operator O is calculated from the link variable 
Uf^{x) instead of U^{x), i.e., {0[U^]) instead of {0[U]). 



We repeat these five steps with different values of the momentum cutoff, and 
observe the dependence of a physical quantity on the momentum cutoff. Then, we 
can determine which energy-momentum component of the gluon field is relevant 
for the physical quantity. Technically, we only have to replace U^{x) instead of 
U^{x), and its operation (mainly the fast Fourier transformation) is numerically 
easy compared to other parts of the simulations. This framework is applicable to 
almost all lattice QCD calculations, and widely useful for QCD phenomena. 

We comment on several points about this framework in the following. 

1. In general, since the gauge transformation is nonlocal in momentum space, the 
energy- momentum component of the gauge field is a gauge-dependent concept. 
The gauge fixing is needed and the result depends on the gauge choice. If we 
do not fix the gauge, the resulting expectation value is zero, even in the case 
of a gauge-invariant operator. 

2. It is possible that the projection in Step 4 contaminates the original condition 
on the momentum cutoff. To evaluate how the projection affects link variables, 
we calculated U^^{x) by repeating Steps 2-4 once again to U^{x), and check 
the overlap between them, iReTr[{?7^(x)}"''f/^^(x)]. The overlap is found to be 
almost unity. For example, the average deviation from unity is about 0.1% at 
Air = 1.5 GeV. Thus, we can expect that the projection does not significantly 
affect the original momentum cutoff. 

3. In general gauges, since the gauge fixing and the momentum cutoff do not 
necessarily commute, it is nontrivial whether these two operations are satis- 
fied simultaneously. In other words, U^{x) can deviate from the gauge fixing 
condition which is originally imposed on U^{x). We have numerically checked 
that, in the case of the Landau gauge, U^{x) almost completely satisfies the 
Landau gauge fixing condition, i.e., Fi,[U^] is maximized. 

4. The Landau gauge is known to suffer from the Gribov copy problem [32] . We 
numerically estimated the systematic uncertainty from the Gribov ambigu- 
ity by random gauge transformation. In the cases of the Wilson loop and 
plaquette, its systematic error is smaller than the statistical error. 
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Chapter 3 

Interquark Potential 



The interquark potential is one of the most fundamental quantities in QCD and 
hadron physics. We here calculate the static quark-antiquark potential. The quark- 
antiquark potential, so-called the Cornell potential [33], is given as 

V{R) =aR-^ + C, (3.1) 

where R is the interquark distance. The physical value of the string tension a is ap- 
proximately 0.89 GeV/fm, and the Coulomb coefficient A is approximately 0.26. The 
constant C depends on the regularization. The quark-antiquark potential includes 
both perturbative and nonperturbative properties. In short range, it is dominated 
by the perturbative one-gluon-exchange Coulomb potential. In long range, it is 
dominated by the nonperturbative linear confinement potential. 

The quark-antiquark potential is extracted from the expectation value of the 
Wilson loop, which is a gauge-invariant path-ordered product of link variables along 
a loop. For statistical improvement, we used the APE smearing method [34j. The 
simulation parameters are listed in Table 13.11 

Table 3.1: Simulation parameters of the gauge configurations. The lattice coupling (3 = 
2Nc/g'^, the lattice volume V, the coordinate-space lattice spacing a, the momentum-space 
lattice spacing Up, and the configuration number A'conf are listed. 





/3 


V 




a [fm] 


ap [GeV] 


-^conf 


Quenched 


5.7 


16^ 


0.19 


0.41 


50 


Quenched 


5.8 


16^ 


0.14 


0.55 


50 


Quenched 


6.0 


16^ 


0.10 


0.77 


50 
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3.1 Interquark potential with the momentum cut- 



We show the quark-antiquark potential with the UV cutoff in Fig. 13.11 When 
the high-momentum gluon is removed by the UV cutoff, the short-range Coulomb 
potential is gradually reduced. This is because the short-range Coulomb potential 
is an UV property. The constant term C is also reduced by the UV cutoff. The 
constant term is mainly given by the lattice regularization for the UV singularity. On 
the other hand, the long-range linear confinement potential is almost unaffected by 
the UV cutoff. At Ayy = 1.5 GeV, the quark-antiquark potential becomes the linear 
potential in the whole region. Note that the quark-antiquark potential would include 
— 7r/(12i?) term, so-called the Liischer term [351 EB]- The Liischer term is similar 
to the perturbative Coulomb potential, but its origin is the nonperturbative string 
ffuctuation. It is nontrivial whether or not the long-range Liischer term remains in 
this calculation. To reveal this, one needs to carefully analyze the functional form 
of the potential in long range with high precision. 

We show the quark-antiquark potential with the IR cutoff in Fig. 13.21 The linear 
confinement potential is strongly affected by the IR cutoff. At Air = 0.1 GeV and 
1.0 GeV, the slope of the linear confinement potential, i.e., the string tension, is 
reduced from the original one. At Air = 1.5 GeV, the linear confinement potential 
completely disappears. The potential becomes the Coulomb-like form. 

From these results, we conclude as follows. The perturbative and nonperturba- 
tive parts of the quark-antiquark potential are decoupled in the momentum space 
of the gluon. The low-momentum gluon induces the confinement potential, and the 
high-momentum gluon induces the Coulomb potential. 

Next, we quantitatively determine which energy- momentum component of the 
gluon field induces the confinement potential. For a quantitative argument, we need 
higher accuracy in momentum space. Since the minimum momentum on lattice is 
the momentum-space lattice spacing a^, the four-momentum length is restricted to 
discrete values as 



We can only take discrete variation on the value of the momentum cutoff. For 
example, the lattice calculations in the range < Air < Op yield the same result. 
This is, so to speak, a discretization error in momentum space. In order to achieve a 
finer resolution in momentum space, we have to calculate with smaller momentum- 
space lattice spacing, i.e., larger coordinate-space lattice volume. 



off 




(3.2) 



3.1. INTERQUARK POTENTIAL WITH THE MOMENTUM CUTOFF 



13 



> 

CD 

CD 
> 



2 - 



6.2 GeV 
4.6 GeV 
3.1 GeV 



Ayy = 1 .5 GeV I — e — i 



/ ^ 

I 

! X 
-S- 













▲ 












▲ 












▲ 




X 


X 




▲ 




M 


X 






A 




M 


B 


B 


B 


B 


M 


B 


B 




e 


e 


e 


B 


e 


e 


e 


1 




1 


0.2 




0.4 




0.6 




0.8 



R[fm] 

Figure 3.1: The static quark-antiquark potential V{K) with the UV cutoff Auv- The 
lattice QCD calculation is performed on 16^ lattice with /3 = 6.0. The broken line is the 
original quark-antiquark potential. 
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Figure 3.2: The static quark-antiquark potential V{R) with the IR cutoff Air. 
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Figure 3.3: The asymptotic string tension (Tasym with the UV cutoff Auv- The vertical 
error bar is the standard statistical error, and the horizontal error bar is the range that 
yields the same result due to the discrete momentum. The original value of the string 
tension is a ~ 0.89 GeV/fm (broken line). 
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Figure 3.4: The asymptotic string tension (Tasym with the IR cutoff Air. 
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For this purpose, we calculate the quark- ant iquark potential on 16^ lattice with 
different three lattice spacings. As listed in Table l3TT| the corresponding momentum- 
space lattice spacings are ~ 0.41 GeV, 0.55 GeV, and 0.77 GeV. We estimate 
the asymptotic string tension cTasym by fitting the quark-antiquark potential with a 
linear function (Tasym-R + const, in 0.3 fm < _R < 0.9 fm. 

We show the asymptotic string tension with the UV cutoff in Fig. I3.3l and with the 
IR cutoff in Fig. 13.41 In these figures, the vertical error bar represents the standard 
statistical error, and the horizontal error bar is not the statistical error but the 
range which yields the same result due to the discrete momentum. In Fig. 13. 3^ the 
asymptotic string tension is almost unaffected in Auv > 1-5 GeV, and its value is 
still its original value a ~ 0.89 GeV/fm. In Auv < 1-5 GeV, the asymptotic string 
tension is drastically reduced. In Fig. 13. 4[ the asymptotic string tension is sensitive 
to the IR cutoff. In Air = 1.5 GeV, the asymptotic string tension is almost zero 
and the confinement potential disappears. 

From these quantitative analyses, we conclude that the relevant energy scale of 
color confinement is below about 1.5 GeV, i.e., color confinement is induced by the 
low- momentum gluon below about 1.5 GeV. 

3.2 Other conditions 

We have calculated with other conditions for the check of consistency: a differ- 
ent gauge and the three-dimensional formalism. Here, we briefly summarize these 
calculations and results. 

Different gauge 

Since our framework is not gauge invariant, it is important to take other gauge 
choices. We calculate with the Coulomb gauge, instead of the Landau gauge. The 
gauge-fixing condition for the Coulomb gauge is to maximize the quantity 

3 

^c[t/]=$^$^ReTrf/,(x). (3.3) 

X j = l 

The fitting results of o"asym are listed in the third column of Table 13.21 The values 
itself are different, but the qualitative behavior is the same as that of the Landau 
gauge. 

Which energy-momentum component induces the confinement potential, in prin- 
ciple, depends on the gauge choice. In the present framework, we can calculate with 
any other gauges in the same way. 
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Tabic 3.2: Asymptotic string tension (Tasym with the IR cutoff Air . The rcsuhs obtained 
under three different conditions arc listed; the main result with the four-dimensional mo- 
mentum cutoff in the Landau gauge, the result in the Coulomb gauge, and the result 
with three-dimensional formalism in the Landau gauge. The lattice QCD calculations are 
performed on 16^ lattice with P = 6.0. 



Air [GeV] q-asyma^(4 - dim.) crasymQ^ (Coulomb) aaayma^ (3-dim.) 

0.051 0.051 0.051 

~ 0.1 0.0469(58) 0.0289(58) 0.0433(51) 

1.0 0.0311(49) 0.0190(59) 0.0198(43) 

1.2 -0.0019(20) 0.0024(25) -0.0034(9) 

1.5 -0.0132(6) -0.0041(10) -0.0092(5) 

2.2 0.0003(12) 0.0058(10) 0.0024(13) 



> 
O 

> 




Figure 3.5: The quark-antiquark potential V{R) with the UV cutoff Auv in the three- 
dimensional formalism. 
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Three-dimensional formalism 

Our framework can be easily extended to the spatial three-dimensional formal- 
ism. We perform the three-dimensional Fourier transformation and consider the 
cutoff by the three- momentum length |p| . The results are listed in the fourth column 
of Table 13. 2[ The behavior is almost the same as the case of the four- dimensional 
momentum cutoff. The quark-antiquark potential with the UV cutoff Auv is shown 
in Fig. 13.51 The confinement potential is not affected in Auv > 1-5 GeV, and 
is affected in Auv < 1-5 GeV. Thus, the confinement potential is induced by the 
low-momentum gluon below about 1.5 GeV also in the three-dimensional formalism. 



3.3 Comparison with an analytical model 

In this section, we analyze the Richardson potential, which is a phenomenological 
model of the interquark potential [2Z|. We introduce the momentum cutoff to the 
Richardson potential and compare it with the lattice QCD result. 

The Richardson potential is constructed so as to reproduce the Coulomb plus lin- 
ear functional structure. It is defined by the one-dressed-gluon-exchange amplitude 
which is proportional to 

Vip') = -Cf'^, (3.4) 

where 

Here, Cp = 4/3 and /3o = (HA'^c — 2A^j)/487r^. A is a parameter in this model. This 
coupling is similar to the standard QCD coupling, except for 1 in the argument of 
logarithm. To obtain the nonrelativistic (instantaneous) potential, we calculate the 
three-dimensional Fourier transformation as 



V{R) = I -^e'R-^V{f 



(27r)3 



■,2^ 1 fiAR) 



(3.6) 



where 



CO g~*^ 
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We set Nf = to compare this potential with the quark- ant iquark potential in 
quenched lattice QCD, and A = 0.48 GeV so that the coefficient of the linear 
potential is equal to the physical string tension a ~ 0.89 GeV/fm. 

We can consider two types of the IR cutoff. One is the simple IR cutoff by the 
three-momentum length on the Fourier transformation, 



V{R) 



JR-p 



Ip1>A:h, (2^)=^ 

The other is the change of the functional form as 



Vip 



{3.t 



V{p 



(3.9) 



IR 



The advantage of the latter way is that we can analytically calculate the momentum 
integral. The result is 



V{R) 



I /A(Ai^) 



(3.10) 



where 



In(A^-l) 
[ln(A2-l)]2+7r2 
1 



ln(l-A2) 



(A>1) 
(1 > A > 0) 



/a(x) 



te 



-tx 



dt- I 



(3.11) 
(3.12) 



and A = Air/A. The symbol V indicates the principal value of the integral, which is 
necessary in the case of A > 1. This functional form becomes the original Richardson 
potential in the limit of Air — > 0, apart from an irrelevant constant. 

These two ways to introduce the IR cutoff yield almost the same results. The 
result of the former way is shown in Fig. 13.61 The irrelevant constant is arbitrarily 
subtracted in the figure. As in the case of lattice QCD, the string tension decreases 
when the IR cutoff is introduced. In Air > 1 GeV, the linear confinement potential 
disappears and the interquark potential becomes a Coulomb-like potential. This 
behavior is consistent with the lattice QCD result. While the Richardson potential 
is only a phenomenological model and its confinement potential is set by hand, it 
well reproduces the AjR-dependence of the interquark potential obtained by lattice 
QCD. 



3.3. COMPARISON WITH AN ANALYTICAL MODEL 
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Chapter 4 
Color Flux Tube 



In a meson, the gluon forms a string-like structure, the so-called color flux tube, 
between the quark and the antiquark. The color flux tube is considered to be 
essential for color confinement. Although the color flux tube is not directly observed 
in experiments, it can be observed in lattice QCD. In this chapter, we apply our 
framework to the analysis of the color flux tube. 

The color flux tube is measured by calculating the action density distribution 
around a static quark-antiquark pair [381 [39] . At the positions of color sources, the 
action density distribution has the peaks which is perturbative singularities of the 
self energies. The color flux tube is formed between the peaks of color sources, but 
its contribution is small compared to perturbative contributions. In Chapter 3, we 
found that the confinement potential originates from the IR gluon below 1.5 GeV 
in the Landau gauge. Thus, if we remove the high- momentum component above 1.5 
GeV by the present framework, we can extract the color flux tube from the action 
density distribution. In other words, we can clearly observe the structure of the 
color flux tube apart from unnecessary perturbative contributions. In this chapter, 
we consider the three-dimensional UV cutoff explained in Section 3.2. 

Table 4.1: Simulation parameters of the gauge configurations. The notation is the same 
as in Table 13. 1[ 
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a [fm] 
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4.1 Action density in vacuum 

Before the action density distribution with a static quark-antiquark pair, we 
calculated the action density without color sources, i.e., in vacuum. The action 
density s{x) is defined as Eq. fl2.7p . The vacuum action density relates to the gluon 
condensate in continuum QCD, which produces the trace anomaly. In the naive 
continuum limit, the action density corresponds to the gluon condensate, however, 
at finite lattice spacing, it is dominated by perturbative corrections [iQl SI]. In 
short, a large part of the vacuum action density on the lattice is perturbatively 
generated. 

The vacuum action density with the UV cutoff is listed in Table W?2\ The 

vacuum action density is translationally invariant, i.e., independent of x. When the 
high-momentum gluon is removed, the vacuum action density drastically decreases. 
For example, at Auv = 1.5 GeV, the vacuum action density is reduced to about 1 % 
of the original value. The remaining small component would lead to nonperturbative 
properties of QCD vacuum. 

Table 4.2: The UV cutoff Auv, the asymptotic string tension ciasym of the quark-antiquark 
potential, and the vacuum action density {s{x)). The statistical error of {s{x)) is omitted 
because it is negligibly small. 



Auv [GeV] [GeV/fm] {s{x)) [a-'] 



No Cut 0.89 14.51 

3.8 0.824(31) 2.57 

2.2 0.801(67) 0.60 

1.5 0.799(18) 0.20 

1.0 0.208(4) 9.5x10^3 



4.2 Action density with a quark-antiquark pair 

The spatial distribution of the action density around a static quark-antiquark 
pair is obtained by measuring s{x) around the Wilson loop at a certain time slice. 
Its expectation value is given by 

. <-W™> - (4.1) 
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where W{R,T) is the value of the Wilson loop with the size RxT. The schematic 
figure is illustrated in Fig. 14.11 The origin of the four-dimensional coordinate is 
placed in the center of the Wilson loop. We measured {s{x))w at the t = plane. 

For statistical improvement, the translationally invariant quantities are averaged 
over, if possible. For example, {p{x)W{R,T)) is given by the convolution sum as 

{pix)WiR, T)) = {^Y1 + ^)^(^' ^' ^))' (4.2) 

s 

where s is the position of the parallel-translated Wilson loop. To enhance the 
ground-state component, the APE smearing method is applied to the spatial link 
variables of the Wilson loop |31]- The action density distribution {s{x))w is indepen- 
dent of the time slice if the ground-state component is suitably dominated. Because 
of the smearing method, the action density distribution is almost independent of T 
in the range of T > 4a in the present calculation. 

We introduce the UV cutoff to the action density distribution as 

(.pA,),, . (£™> _ ,,3, 

The arguments, such abbreviated for simplicity. The physical interpretation 

is the spatial distribution of the low-momentum gluon around a physical quark- 
antiquark pair. As another choice, one can introduce the UV cutoff not only for 
s[U] but also for VI^[?7], and its result is expected to be qualitatively similar to 
Eq. (lOj). 

In Fig. 14. 2[ we display the action density distribution {s{x))w with the UV 
cutoff Auv = 1.5 GeV, 2.2 GeV, and 3.8 GeV. The interquark distance between 
the quark and the antiquark is R = 0.6 fm. The overall sign of {s{x))\y is flipped 
in the flgure, which is only a matter of deflnition. Because the absolute value of 
{s{x))w at Auv = 1-5 GeV and 2.2 GeV is small, these data in Fig. 14.21 are enlarged 
by a factor of flve compared to the other ones. The action density distribution in 
original lattice QCD without the momentum cutoff is also displayed ("No Cut" in 
the upper left); however, its statistical error is relatively large. In original lattice 
QCD, the action density is strongly enhanced at the positions of the quark and 
the antiquark. In contrast, the color flux tube is difficult to observe because of 
such singular peaks and the large statistical fluctuation. Both the singular peaks 
and the large fluctuation originate from the perturbative part of the action density. 
When the high-momentum gluon above 3.8 GeV is removed (the upper right), these 
perturbative contributions are drastically suppressed, and the flux-tube structure 
connecting the quark and the antiquark becomes clear. At Auv = 1-5 GeV and 
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R 

< > 




Figure 4.1: The Wilson loop W{R,T) and the three-dimensional plane where the action 
density {s{x))w is measured. The origin of the coordinate is placed in the center of the 
Wilson loop. 




Figure 4.2: The action density distribution {s{x))w around a static quark-antiquark pair. 
The separation between the quark and antiquark is i? = 0.6 fm. "x5"s mean that the 
lower figures are five times enlarged in the vertical direction compared to the upper figures. 
The "No Cut" is the original lattice QCD resuh. 
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2.2 GeV (the lower right and the lower left, respectively), the two peaks seem to 
disappear, and the action density is distributed around the origin. 

Apart from the vacuum contribution, which is translationally invariant, the ac- 
tion density at Auv = 1-5 GeV is broadly distributed around the midpoint between 
the quark and the antiquark. In the calculation of the quark-antiquark potential in 
Chapter 3, this low-momentum gluon leads to the linear confinement potential over 
the entire range of R. Therefore, this action density distribution corresponds to the 
color flux tube inducing the confinement potential. 

In Fig. 14.31 we plot the action density distribution {s{x))w along the x-axis, 
i.e., in the longitudinal direction of the quark-antiquark separation. In original 
lattice QCD without the UV cutoff ("No Cut"), the action density has two self- 
energy peaks at the positions of the quark and the antiquark, x = 0.3 fm and 
X = —0.3 fm, respectively. When the UV cutoff is introduced, these self-energy 
peaks are drastically suppressed. Moreover, the absolute values of the action density 
and the statistical fluctuation become small, as the vacuum action density. The 
results at Auv = 1-5 GeV and 2.2 GeV are also shown in Fig. 14. 4[ The action 
density distribution has a maximum at the origin, and the self-energy peaks seem 
to disappear. Although the self-energy peaks would also include a nonperturbative 
contribution, it is too small to distinguish from the flux tube. The endpoints of 
the flux tube are not sharp, spreading outside the positions of the quark and the 
antiquark. In the case of R = 0.6 fm, the longitudinal shape of the flux tube 
resembles a broad mountain rather than a plateau. 

We also show the energy density distribution {e{x))w along the x-axis in Fig. 14.41 
The energy density e{x) is defined identically as the action density by changing 
the relative sign between the spatial plaquettes and the temporal plaquettes. The 
absolute value of the energy density is smaller than that of the action density due to 
the cancellation between the chromoelectric contribution and the chromomagnetic 
contribution. Apart from the absolute value, the overall shape of the energy density 
distribution is similar to that of the action density distribution. 

To estimate the width of the flux tube, we fit the action density distribution 
along the y-axis to a the Gaussian form, sq exp(— y^/5^). It is seen that the Gaussian 
form can well reproduce the lattice data, as shown in Fig. 14.51 The best-fit width 
parameter 6 is 0.31±0.01 fm at Auv = 2.2 GeV, and 0.35±0.01 fm at Auv = 1-5 
GeV. These values are comparable to the fiux-tube width of earlier works in the 
standard lattice QCD [121 SSI SI SS] . 

Next, we analyze how the fiux-tube shape depends on the interquark distance R. 
As shown in Fig. 14.21 the longitudinal length and the transverse width are almost the 
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Figure 4.3: The action density distribution {s{x))\\r along the x-axis. The interquark 
distance R is 0.6 fm. The quark and the antiquark are located on x = 0.3 fm and 
X = —0.3 fm, respectively. The "No Cut" is the original lattice QCD result. 
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Figure 4.4: The action density distribution {s{x))w and the energy density distribution 
{e{x))w along the x-axis. The interquark distance R is 0.6 fm. The quark and the 
antiquark are located on a; = 0.3 fm and x = —0.3 fm, respectively. 
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Action Density Ayv=2. 2 GeV i — • — i 
Energy Density Auv=2.2 GeV i — e — i 
Action Density Ayy=1 .5 GeV i — ■ — i 
Energy Density Auv=1 -5 GeV i — b — i 




y[fm] 

Figure 4.5: The action density distribution {s{x))w and the energy density distribution 
{£{x))w along the y-axis. The interquark distance R is 0.6 fm. The sohd hues are the 
results of fitting to a the Gaussian function sq exp(— y^/(5^). 




Figure 4.6: The action density distribution {s{x))y[r for the interquark distance R = 1.0 
fm. This calculation is performed on a 32x16^ lattice. 
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same at R = 0.6 fm and Auv = 1-5 GeV. The overall shape seems to be isotropic 
and far from a "string" or "tube". This is because the transverse width of the 
flux tube is fairly large. To approach the distribution to the tube-like shape, the 
interquark distance must be enlarged [l3l HH HS]. Here, we use a 32(x-axis)xl6^ 
lattice instead of a 16^ lattice and extend the interquark distance to R = 1.0 fm. 
As shown in Fig. 14. 6^ the distribution is stretched in the longitudinal direction, and 
the longitudinal length becomes larger than the transverse width. Compared to the 
result of R = 0.6 fm, the flux tube at i? = 1.0 fm approaches a broad tube-like 
shape. 
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Chapter 5 
Meson Mass 



In this chapter, we calculate the mass of a pion (pseudoscalar meson) and a p- 
meson (vector meson) in quenched lattice QCD. The physical values are ~ 140 
MeV and rrip ^ 770 MeV. The ground-state mass of a meson is extracted from the 
large-time behavior of the corresponding meson correlator as 

{Ml{x,t)Mr{x,0)) = Coe-'"^* + ■ ■ ■ . (5.1) 

The interpolation operator Mr(x, t) is, for example, 9755 for a pion and q'jjq for a 
p-meson. The dots contain the excited-state components, which rapidly die out in 
large t. 

We used the clover fermion and the staggered fermion for the quark propagator. 
The clover fermion is the 0(a)-improved Wilson fermion, and the quark mass is 
described in terms of the mean- field-improved hopping parameter k [IHl HZ]. The 
setups of the gauge configurations are shown in Table 15.11 The clover fermion and 
the staggered fermion are used on 16^ and 16^ x 32 lattices, respectively. 

Table 5.1: Simulation parameters of the gauge configurations. The notation is the same 
as in Table [XH 
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5.1 Meson mass with the momentum cutoff 

In Fig. 15. H we show the meson masses with the UV cutoff with the clover fermion. 
The data at Auv — 12.5 GeV are the original lattice results. As the high-momentum 
components are removed, both the pion mass m^r and the p-meson mass rrip gradually 
decreases. This is because the UV gluon "dresses" the quarks in the mesons. 

In Fig. 15.21 we show the meson masses with the IR cutoff with the clover fermion. 
As in the case of the UV cutoff, the meson masses are reduced by the IR cutoff. 
In addition, the pion mass and the p-meson mass degenerate in Air > 1.5 GeV. 
In Chapter 3, we found that the confinement potential disappears in this region. 
Thus, this degeneracy suggests that the quark and antiquark become unbound or, if 
possible, very narrowly bound. In such a state, the character of each meson is lost, 
and the meson mass is equal to the sum of the quark and antiquark masses. This 
state is called as "quasi- free" . 

Next, we consider the dependence of the meson masses on the quark mass with 
the staggered fermion. For chiral extrapolation, the squared meson masses are fitted 
with quadratic functions as 



where C2, Ci, and Cq are fitting parameters. We show the results of chiral extrapola- 
tion in Fig. 15.31 and Fig. 15.41 When the IR cutoff is introduced, the p-meson mass 
uniformly decreases but the linear extrapolation form is unchanged. In contrast, 
the functional form of the pion mass changes in Air > 1.5 GeV, where the pion 
mass degenerates the p-meson mass. The extrapolation function becomes a linear 
function. 

The pion mass is related to the chiral condensate through the Gell-Mann-Oaks- 
Renner relation, 



In principle, we can examine spontaneous chiral symmetry breaking through the 
behavior of the pion mass. However, especially when the quarks are deconfined, 
we have to take care of the finite- volume artifact. We directly analyze spontaneous 
chiral symmetry breaking by calculating the chiral condensate in the next chapter. 



(m^a)^ or {nipof' = C2{maf' + Ci(ma) + Cq, 





m{qq). 



(5.3) 
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Figure 5.1: The Auv-dependence of pion mass rriTr and p- meson mass nip with the clover 
fermion. The data at Auv — 12.5 GeV is the original lattice result. 
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Figure 5.2: The AiR-dependence of meson masses with the clover fermion. The data at 
Air = GeV is the original lattice result. 
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Figure 5.3: The rho mass rup plotted against the bare quark mass m with the staggered 
fermion. The "No Cut" (broken hne) is the original lattice QCD result. 
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Figure 5.4: The pion mass m^r plotted against the bare quark mass m with the staggered 
fermion. 



32 



CHAPTER 6. CHIRAL CONDENSATE 



Chapter 6 

Chiral Condensate 



In this chapter, we analyze the chiral condensate. The chiral condensate is an 
order parameter of spontaneous chiral symmetry breaking in the chiral limit. It is 
nonzero in the symmetry-broken phase and zero in the symmetry-restored phase. 
When the quark mass is finite, the chiral condensate includes the effect of explicit 
breaking by the quark mass as well as the dynamical effect. To extract the chiral 
limit in lattice QCD, one calculates with several quark masses and extrapolates to 
the chiral limit. We denote the absolute value of the chiral condensate as 

J: = -a^qq) = ahrSq, (6.1) 

where Sg is the quark propagator. 

We compute the chiral condensate with the staggered fermion in quenched QCD 
and in full QCD. The full QCD configuration includes the two-flavor staggered 
fermion. In full QCD, we use a single mass for the valence and sea quarks, ma = 
ruscaa = 0.01, of which the corresponding pion mass is about 500 MeV. The simula- 
tion parameters of gauge configurations are summarized in Table 16. 1[ 

Table 6.1: Simulation parameters of the gauge configurations. The notation is the same 
as in Table [STTl The full QCD configuration is provided in the NERSC archive [48j. rrisea. 
is the bare quark mass of the staggered fermion. 
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6.1 Chiral condensate with the UV cutoff 

First, we show the chiral condensate with the UV cutoff in Fig. 16.11 Since there 
is no significant difference between the quenched and full QCD results, we plot only 
the full QCD result. The right-side point at Ayv — 12.5 GeV is the result of original 
lattice QCD without the momentum cutoff. 

Although spontaneous chiral symmetry breaking is expected to be caused by 
nonperturbative gluons, the chiral condensate is drastically changed by the UV cut- 
off. However, as shown below, this is mainly because the chiral condensate is a 
renormalization-group variant and UV-diverging quantity. It is dressed by pertur- 
bative gluons and its value strongly depends on the UV regularization. In standard 
lattice QCD, the perturbative contribution is several orders of magnitude larger than 
the nonperturbative core of the condensate, as shown in Section 4.1. 

To estimate the effect of renormalization, we calculate a renormalization factor, 
so-called a Z-factor, nonperturbatively ^\ \5U\. The renormalization factor Zo{k) 
is determined from the amputated Green function of the quark bilinear operator O. 
The renormalization condition is imposed as 

Zo{k)Z;\k)To{k) = 1, (6.2) 

where 

^o{k) ^ Y^tT[S;\k)Go{k)S;\k)Pl] (6.3) 

Go{k) ^ {q{k)Om) (6.4) 
S,{k) ^ {q{k)q{k)). (6.5) 

Pq is the appropriate projection operator. The wave-function renormalization factor 
Zy^(fc) of the quark field is obtained from the conserved vector current, i.e., Zv{k) = 
1. Note that k is the momentum of the quark field, not the momentum of the gluon 
field. 

We calculate the renormalization factor Zs{k) of the scalar operator O = qq, and 
plot the renormalized chiral condensate Zs{5 GeV) x S in Fig. I6.1[ The renormalized 
chiral condensate is almost independent of the UV cutoff. As the UV gluon is 
removed by the UV cutoff, the bare chiral condensate approaches the renormalized 
one. This means that the drastic change by the UV cutoff is well explained in terms 
of renormalization. 
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Figure 6.1: The chiral condensate S = —a^{qq) with the UV cutoff Auv The quark mass 
is ma = 0.01. The "renormalized" chiral condensate is multiphed by the renormahzation 
factor Zs- 

6.2 Chiral condensate with the IR cutoff 

Second, we analyze the chiral condensate with the IR cutoff. We show the full 
QCD result in Fig. 16.21 and the quenched QCD result in Fig. 16.31 The quark mass is 
ma = 0.01 in both calculations. In the case of the IR cutoff, the chiral condensate 
does not show the drastic change. Thus, we expect the physical contribution to 
spontaneous chiral symmetry breaking instead of an artifact of renormahzation. 

When the IR gluon is removed, the effective quark mass would be reduced espe- 
cially at a large distance. Thus, we must pay attention to the finite-volume effect 
in Air > 0, even though our lattice volume is large enough at Air = 0. We estimate 
the finite-volume effect by changing boundary conditions of the quark propagator 
|51] . In Fig. 16.21 and Fig. 16. 3[ "PBC" and "APBC" mean periodic and antiperiodic 
boundary conditions, respectively. Since the result is independent of the boundary 
conditions if the lattice volume is large enough, the difference between these data 
should be understood as the finite- volume effect. As seen from Fig. 16. 2[ the 16^ x 32 
lattice of full QCD suffers from the finite- volume effect in Air > 1.0 GeV. From 
Fig. 16. 3[ the finite- volume effect is fairly small for the 32^ lattice of quenched QCD, 
although it gradually grows in Air > 1.5 GeV. 

Both in Fig. 16.21 and Fig. 16. 3[ the chiral condensate suddenly gets small around 
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Figure 6.2: The full QCD result of the chiral condensate S = —a^{qq) with the IR cutoff 
Air. The quark mass is ma = 0.01. PBC and APBC mean periodic and antiperiodic 
boundary conditions, respectively. 
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Figure 6.3: The quenched QCD result of the chiral condensate with the IR cutoff. 



36 



CHAPTER 6. CHIRAL CONDENSATE 



Mr = 0. This jump around Air = is caused by cutting only the zero-momentum 
hnk variable f/^(0). Despite the change at a single point = 0, the chiral con- 
densate is about 40% reduced. Such a large change is not observed in removing 
other low-momentum components. Therefore, the zero-momentum gluon is special 
and it possesses a major contribution to the chiral condensate. Note that "zero 
momentum" on momentum-space lattice corresponds to the deep-IR region which 
is roughly < in the continuum. 

The zero-momentum link variable corresponds to a spatially-uniform gauge back- 
ground. In general, the non-Abelian gauge field could have a nontrivial effect even in 
spatially-uniform case, unlike the Abelian gauge field. Our result actually suggests 
that the zero-momentum gauge field contributes to the chiral condensate. Note, 
however, that it is nontrivial whether spontaneous chiral symmetry breaking occurs 
only by the spatially-uniform gauge background. 

In large Air, since the lattice volume of full QCD is not large enough, we analyze 
the quenched QCD result in Fig. 16.31 When the "nonzero-momentum" gluon of 



a/p^ > ttp is removed by the IR cutoff, the chiral condensate gradually decreases. 

Thus, not only the "zero-momentum" gluon but also the "nonzero-momentum" 
gluon contributes to the chiral condensate. The chiral condensate continues to 
decrease even in Air > 1.5 GeV. Although it is difficult to perform an accurate 
analysis in large Air because of the finite- volume effect, we can see that the chiral 
condensate is affected by the gluon in the intermediate-momentum region. 

Next, we consider the chiral extrapolation of the chiral condensate. When the 
bare quark mass m is small, the chiral condensate is expanded as a function of m, 
as 



where S'(m) = dTi{m)/dma. S(0) represents spontaneous chiral symmetry breaking 
in the chiral limit. We fit the quenched QCD result by the linear extrapolation 
function E(0) +maE'(0). The fitting result is shown in Fig. 16.41 and Table [621 Note 
that the data of "Air ~ 0.1 GeV" corresponds to the smallest IR cutoff, which cuts 
only the zero- momentum link variable, and so the value "0.1 GeV" itself is not so 
meaningful. 

As stated above, when the zero-momentum gluon field is removed, the chiral 
condensate is largely changed. S(0) is about 40% reduced and S'(0) is about 30% 
reduced. As for the nonzero-momentum gluon, the extrapolating line moves down 
parallel by the IR cutoff. S(0) is gradually reduced and S'(0) is almost unchanged. 
This indicates that the nonzero-momentum gluon has small but finite contribution 
to spontaneous chiral symmetry breaking. 




S(m) = S(0)+maS'(0) + -- - , 



(6.6) 
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Figure 6.4: The chiral extrapolation of the chiral condensate S = —a'^{qq). Air ~ 0.1 
GeV corresponds to the cutoff for the zero-momentum hnk variable. 

In Fig. I6.4[ our result suggests another interesting possibility. At least within the 
present numerical accuracy, the chiral condensate in the chiral limit remains finite 
at Air = 1.5 GeV, which is the relevant energy scale of color confinement. If this 
is true, this means that the relevant energy scale of spontaneous chiral symmetry 
breaking is larger than that of color confinement at zero temperature. Unfortunately, 
however, we cannot make a decisive statement because of systematic error of the 
chiral extrapolation. For more conclusive answer, we need the full QCD calculation 
very close to the chiral limit, while the finite- volume effect is severely crucial in large 
Air and small m. 

The Banks-Casher relation states that the chiral condensate is related to the 
spectral density p(A) of the Dirac operator as 

E = 7rp(0), (6.7) 

in the chiral limit [13] . The spectral density of the Dirac operator is given in infinite 
volume as 

p(A)= lim i^5(A-A,), (6.8) 

k 
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Table 6.2: The fitting result of the chiral extrapolation in Fig. 16.41 The extrapolation 
function is S(0) + maS'(O). 

Air [GeV] S(0) S^O) 

0.00639(81) 1.269(53) 

~ 0.1 0.00380(27) 0.933(16) 

1.5 0.00200(7) 0.948(4) 

1.8 0.00155(2) 0.929(1) 



and the eigenvalue of the Dirac operator is iXk- The zero mode of the quark field is 
directly related to spontaneous chiral symmetry breaking from this relation. In con- 
trast, the gluon field is nontrivially related to spontaneous chiral symmetry break- 
ing. Our result presents the connection between the energy-momentum component 
of gluons and the zero mode of quarks, i.e., the energy scale of gluons and quarks. 
We directly analyze the Dirac spectrum in the next chapter. 



39 



Chapter 7 
Dirac spectrum 

In this chapter, we analyze the eigenvalue spectrum of the Dirac operator. Al- 
though the Dirac equation is written in a simple form, the Dirac spectrum tells 
us nontrivial properties of QCD. In particular, zero modes of the Dirac operator 
have special roles on nonperturbative properties induced by the gluon field, e.g., 
spontaneous chiral symmetry breaking and topological charge. Moreover, the low- 
lying Dirac spectrum shows the universal behavior of disordered systems, which is 
described by chiral random matrix theory. It is important to clarify how the gluon 
field induces these properties of the Dirac spectrum. 

We computed the low-lying eigenvalues of the massless staggered Dirac operator 
in quenched lattice QCD. The parameters of the gauge configurations are listed 
in Table I7.1[ We here adopted the 1-loop improved Symanzik gauge action with 
tadpole improvement [52j. The tadpole coefficient is defined as the fourth-root of 
the average plaquette value [17|. We generated four kinds of gauge configurations. 
The 12^, 16^, and 20^ lattices are used to analyze the eigenvalues and the chiralities, 
and the 8^ lattice is used to discuss the connection with chiral random matrix theory. 

Table 7.1: Simulation parameters of the gauge configurations. The notation is the same 
as in Table 13. 1[ 





/3 


V 


[A 


a [fm] 


ttp [GeV] 


A'conf 


Quenched (improved) 


9.0 


20^ 


0.07 


0.89 


50 


Quenched (improved) 


8.6 


16^ 


0.08 


0.97 


50 


Quenched (improved) 


8.3 


12^ 


0.11 


0.94 


50 


Quenched (improved) 


7.9 


8 


4 


0.16 


0.97 


5000 
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7.1 Eigenvalue and chirality 



The j-th eigenmode of the Dirac operator is given by 



Dijjj = iXjipj- 



(7.1) 



Because the massless staggered Dirac operator is anti-Hermitian, the eigenvalues 
always appear in pairs, ±iXj. We only have to calculate the positive eigenvalues Xj > 
0. It should be understood that the negative eigenvalues with the same magnitude 
exist. 

As explained in Section 2.1, the lattice staggered fermion breaks the flavor sym- 
metry, i.e., the full SU(A^/) chiral symmetry. As a result, the unimproved staggered 
Dirac operator fails to reproduce nontrivial topological sectors [531 EH ES]. This 
problem is approximately settled by improving the lattice gauge action and the lat- 
tice fermion action [561 EH EH EH]. For the improved gauge action, we adopted the 
1-loop improved Symanzik gauge action with tadpole improvement [S2]- For the 
improved staggered Dirac operator, we adopted the FATTxASQ operator The 
FATTxASQ operator is constructed from the FAT7 link variable and the ASQ oper- 
ator [61]. The spatial boundary conditions are periodic, and the temporal boundary 
condition is antiperiodic. 

In the top of Fig. 17. H we show the low-lying 10 positive eigenvalues of one 
typical gauge configuration of /3 = 8.6. This improved operator strongly suppresses 
the flavor-symmetry breaking lattice artifact. The four-fold degeneracy, i.e., the 
four-flavor symmetry, is satisfied with high precision. In addition, zero eigenvalues 
are approximately reproduced. Such zero modes are called as "would-be" or "near" 
zero modes. Because the negative eigenvalues with the same magnitudes exist, the 
topological charge of this gauge configuration is Q = 1. 

In Fig. I7.1[ we also show the Dirac spectrum of the same gauge configuration with 
Auv = 3 GeV and Ajr = 3 GeV. For Auv = 3 GeV, while all the eigenvalues change, 
the four-fold degeneracy is still satisfied. The would-be zero modes disappear. This 
means that the topological charge is broken. For Air = 3 GeV, the low-lying 10 
eigenvalues become close to the eigenvalue of free fermions, which is about 0.39 for 
this case. Also in this case, the topological charge becomes zero. 

Next, we discuss the correlation between the eigenvalue and the chirality. The 
chirality is defined as 



For nonzero eigenvalues, their chiralities are zero. On the other hand, for zero 
eigenvalues, their chiralities are +1 or —1. In the top of Fig. 17.21 we show the 



(7.2) 
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Figure 7.1: The eigenvalue spectrum of the FAT7xASQ Dirac operator. Only low- lying 
10 positive eigenvalues of one typical gauge configuration are shown. 



scattering plot of all the low-lying 30 eigenvalues of 50 gauge configurations of /9 = 
8.6. There exists a clear separation between the nonzero modes, which is \xj\ — 0, 
and the would-be zero modes, which is \xj\ — 0.8. The reason why the chirality 
deviates from ±1 is that it requires a renormalization |62] . 

In the middle of Fig. \7.2\ we show the scattering plot with the UV cutoff Auv = 
6, 4, and 3 GeV. As the high-momentum component is removed, the eigenvalue 
spectrum is gradually changed. For Auv = 3 GeV, all eigenmodes are \xj\ — 0. 
Therefore, all gauge configurations fall into the trivial topological sector. In the 
bottom of Fig. 17.21 we show the scattering plot with the IR cutoff Air = 1, 2, 
and 3 GeV. When the low-momentum component is removed, many of the low- 
lying eigenvalues approach the free eigenvalue, which is \j ~ 0.39 and \xj\ — 0. 
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Figure 7.2: The scattering plot of the eigenvalue Xj and the absolute value of the chirality 
Xj- The low-lying 30 eigenvalues of 50 gauge configurations are plotted. The calculation 
is done on 16^ lattice with /? = 8.6. 
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Compared to the case of the UV cutoff, the changes of the eigenvalues are sensitive 
to the IR cutoff. This resuh is consistent with our expectation that the low-lying 
Dirac eigenmode interacts with the low-momentum gluon. For Air = 3 GeV, we 
cannot identify would-be zero modes since all eigenvalues are almost degenerate. 
From these results, except for the very-high-momentum components, the relatively 
broad region of momentum components of the gluon field is related to the low-lying 
Dirac spectrum. 

When the momentum cutoff is imposed, the clear separation between would-be 
zero modes and nonzero modes is lost. This is because the staggered Dirac operator 
does not have exact zero modes. This situation will be improved by the overlap 
operator, which has exact zero modes [29t [63]. 



7.2 Topological charge 

The numbers of zero modes determine the topological charge Q of the background 
gauge field through the Atiyah-Singer index theorem [61] . The index theorem states 

ur -nL = NjQ, (7.3) 

where nn and ul are the numbers of right-handed and left-handed zero modes, 
respectively [5^ . 

From the number of would-be zero modes, we extracted the topological charge Q. 
We set the criterion for would-be zero modes as \xj\ > 0.6, and for nonzero modes as 
\Xj\ < 0.6. This criterion clearly separates would-be zero modes in the original case, 
as shown in the top of Fig. 17. 2[ For the cases in which the clear separation between 
would-be zero modes and nonzero modes is lost, this criterion could be subtle. The 
results are summarized in Table 17. 2[ 

For a consistency check, we also extracted the topological charge in another way. 
We directly calculated the topological charge by 

Q = ^ f wTr[F^.(x)F.p(x)]. (7.4) 

fiuap 

Its discretized form is 

1 

Q = -^j^E E WTr[f/;..(x)f/.p(a:)], (7.5) 

where Ufj^^^x) is the plaquette. We measured this topological charge after the APE 
smearing steps [31] 155] . The smearing coefficient and the number of the smearing 
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step are set to 7.0 and 200, respectively. The results are shown in the parentheses in 
Table 17.21 The obtained topological charge is roughly consistent with the estimate 
by the would-be zero modes, especially in the original lattice QCD, as expected. 
These two estimates deviate from each other in some cases with the momentum 
cutoff, e.g., Auv = 4 GeV. This is partly because it is difficult to identify the would- 
be zero modes in these cases, and partly because the APE smearing steps change 
the local topological structure. 

We can see that the topological charge is destroyed by removing the broad mo- 
mentum region of the gauge field. This behavior is seen in the both cases of the UV 
and IR cutoffs. Thus, there is no typical energy scale of the gauge field inducing 
the topological charge. This is different from the case of color confinement, which is 
induced only by the narrow low-momentum component of the gauge field, as shown 
in Chapter 3. 

In Table 17. 3^ we show the numerical results calculated in the Coulomb gauge, 
instead of the Landau gauge. The results are almost consistent with those of the 
Landau gauge. 



7.3 Connection with chiral random matrix theory 

Another interesting property of the low-lying Dirac spectrum is that it is suc- 
cessfully described by chiral random matrix theory. Chiral random matrix theory 
is an effective theory which describes the universal behavior of disordered systems 
[65| [66| [67] . Chiral random matrix theory suggests that the low-lying Dirac spec- 
trum depends only on the global symmetries of the system, and not on the details of 
the microscopic structure. Once the validity of chiral random matrix theory is guar- 
anteed, chiral random matrix theory enables us to predict the low-lying eigenvalue 
distribution by a simple and universal function. 

We compared the lattice data with a prediction of chiral random matrix theory, 
using the gauge configurations of (3 = 7.9. We analyzed the Q = sector with the 
criterion for nonzero modes as \xj\ < 0.1. 

We calculated the probability distribution of the lowest nonzero eigenvalue Amin- 
In chiral random matrix theory, the probability distribution of the lowest eigenvalue 
is given as 

P(Amin) = (7.6) 

^min = AminS^ {"^ ■'^) 

for the SU(3) quenched gauge field with Q = [68] • If one directly calculates the 
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Table 7.2: Topological charge Q estimated by would-be zero modes. The estimate by 
Eq. ()7.5p is also shown in the parenthesis. The results of the 50 gauge configurations with 
/? = 8.6 are shown. 



V [a^] 


Data 


\Q\ 


= 


IQI = 1 


\Q\ 


= 2 


\Q\ 


= 3 




No Cut 


28 


(29) 


18 (17) 


4 


(4) 





(0) 




Auv = 6 GeV 


35 


(27) 


12 (17) 


3 


(6) 





(0) 




Auv = 4 GeV 


50 


(44) 


0(6) 





(0) 





(0) 




Auv = 3 GeV 


50 


(50) 


0(0) 





(0) 





(0) 




Air = 1 GeV 


29 


(27) 


17 (19) 


4 


(4) 





(0) 




Air = 2 GeV 


32 


(43) 


13 (7) 


5 


(0) 





(0) 




Air = 3 GeV 


49 


(50) 


1 (0) 





(0) 





(0) 




No Cut 


28 


(28) 


22 (22) 





(0) 





(0) 




Auv = 4 GeV 


50 


(50) 


0(0) 





(0) 





(0) 




Auv = 3 GeV 


50 


(50) 


0(0) 





(0) 





(0) 


12^ 


No Cut 


28 


(28) 


20 (20) 


2 


(2) 





(0) 




Auv = 4 GeV 


49 


(44) 


1(6) 





(0) 





(0) 




Auv = 3 GeV 


50 


(50) 


0(0) 





(0) 





(0) 



Table 7.3: Topological charge Q in the Coulomb gauge. 



V [g^] Data \Q\ = \Q\ = 1 \Q\ = 2 \Q\ = 3 

16^ No Cut 28 (29) 18 (17) 4 (4) (0) 

Air = 1 GeV 29 (27) 17 (20) 4 (3) (0) 

Air = 2 GeV 35 (31) 13 (14) 1 (5) 1 (0) 

Air = 3 GeV 49 (50) 1 (0) (0) (0) 



chiral condensate in infinite volume, this function is parameter free. In this study, 
we treated E as a fitting parameter. The lattice data and the best-fit result of 
Eq. (17. 6p are shown in Fig. 17.31 In the original lattice QCD, the lattice data is well 
reproduced by chiral random matrix theory. The best-fit value is S ~ 0.0022 x a~^. 

For Auv = 4 GeV, the qualitative behavior does not change. The lattice data is 
still reproduced by Eq. fl7.6p . Thus, the high- momentum gluon is irrelevant for the 
universality of the lowest eigenvalue distribution. The best-fit value of E is changed 
by the UV cutoff. This is because the chiral condensate is renormalization-group 
variant and UV divergent. By definition, the UV divergent quantity depends on the 
UV cutoff, although this might be irrelevant for phenomenology. 

For Air = 2 GeV, the lowest eigenvalue distribution becomes to concentrate 
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on the vicinity of the free eigenvalue. The lattice data drastically deviates from 
Eq. (17. 6p . Therefore, the low-momentum gluon is crucial for the validity of chiral 
random matrix theory. In other words, the low-momentum component induces the 
strong-interacting and disordered nature of the gauge field. 

In the calculation of the chiral condensate in Chapter 6, the zero- momentum com- 
ponent has a large contribution to the chiral condensate. To check it, we calculated 
the Dirac spectrum with the small IR cutoff Air ~ 0.1 GeV. For Air ~ 0.1 GeV, the 
lattice data can be fitted by Eq. (17. 6p with the best-fit value S ~ 0.0019 x a'^. The 
contribution of the zero-momentum component is indeed large, about 15% of the 
total. This is qualitatively consistent with the result in Chapter 6. Note, however, 
that this is quantitatively different because the lattice action and the lattice spacing 
are different. 
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Figure 7.3: The probability distribution of the lowest eigenvalue Amin in the Q = sector. 
The histograms are the lattice data of 8^ lattice with /3 = 7.9. The solid curves for No 
Cut and Auv = 4 GeV are the best-fit functions of Eq. ()7.6p from chiral random matrix 
theory. 
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Chapter 8 
Summary 

We have introduced the lattice framework to analyze the relation between a 
QCD phenomenon and the energy-momentum component of the gluon field. In 
this framework, using the Fourier transformation of the link variable, we construct 
the momentum-space link variable, remove some energy-momentum components, 
and reconstruct the coordinate-space link variable. This framework is applicable to 
many lattice QCD calculations. We have applied this framework to the analysis of 
color confinement and spontaneous chiral symmetry breaking. 

As for color confinement, we have analyzed the quark- ant iquark potential, the 
color flux tube, and meson masses. From the quantitative analysis of the quark- 
antiquark potential, we have found that the confinement potential is induced by 
the low-momentum component below 1.5 GeV in the Landau gauge. The low- 
momentum component below 1.5 GeV produces the linear interquark potential, 
and the high- momentum component above 1.5 GeV produces the Coulomb-like in- 
terquark potential without confinement. As one application of this result, by re- 
stricting the gluon field to the high- momentum component above 1.5 GeV, we can 
extract the color flux tube from the action density distribution. In addition, when 
this low-momentum component is removed, the two quarks in mesons are not con- 
fined and become quasi-free. 

As for spontaneous chiral symmetry breaking, we have analyzed the chiral con- 
densate and the Dirac spectrum. The chiral condensate is induced by the broad 
low- momentum region, even above 1.5 GeV. This behavior is also observed in the 
analysis of the Dirac spectrum. This dependence on the momentum cutoff is different 
from that of color confinement. The present result suggests that spontaneous chiral 
symmetry breaking is induced by higher energy-momentum component of gluons, 
compared to color confinement. This behavior is consistent with phenomenological 
models [TH [T5] . For more conclusive statement, we have to take into account many 
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systematics, especially, the sea quark effects in the case of realistic QCD. 

We have also analyzed topological charge in the two ways: the fermionic way 
through the Dirac zero mode and the field theoretical way through the gauge field 
distribution. The momentum cutoff gradually affects the Dirac spectrum and the 
zero modes, and then changes the topological charge. This is similar to the behavior 
of the chiral condensate. 

In this framework, we need gauge fixing, and the obtained result depends on the 
gauge choice. Nevertheless, this analysis is important in the following theoretical 
viewpoints. 

First, the relevant energy scale, i.e., the energy-momentum component inducing 
a QCD phenomenon, would be useful for developing effective theories. An effective 
theory has some cutoff to simplify physics or to reduce degrees of freedom. Although 
the value of the cutoff should be based on some physical reasoning, its microscopic 
derivation is difficult in many cases. The relevant energy scale obtained by lattice 
QCD would provides a physical reasoning for the cutoff value. For example, by 
setting the UV cutoff to be 1.5 GeV in momentum integral, we can safely pick up 
the contribution to color confinement. Further, the relevant energy scale determines 
degrees of freedom which appear in the effective field theory |69] . 

Second, this analysis reveals the relation between different QCD phenomena 
in terms of the energy-momentum component of the gauge field. We studied the 
connection between color confinement and spontaneous chiral symmetry breaking. 
Such connection would be important for understanding their underlying mecha- 
nisms. Also, based on the connection between the confinement potential and the 
color flux tube, we succeeded in extracting the color flux tube from the gluon dis- 
tribution. 

Third, this framework has a technical advantage in numerical simulation. As 
shown in Chapter 4, the statistical fluctuation of the action density is strongly sup- 
pressed by removing the UV components. This is because the statistical fluctuation 
is mainly given by the high-frequency mode, i.e., the UV component of the link vari- 
able. The statistical error is one major difficulty in lattice QCD. This framework 
can reduce the statistical error without changing the target IR physics [70] . 
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